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We calculate the electron-energy loss spectrum and the optical conductivity for NaV2 0s using 
the standard Lanczos algorithm. The vanadium ions in NaV2 0s form a system of coupled ladders 
which can be described by a quarter-filled extended Hubbard model. Since this system has a large 
unit cell, one has to be very careful to avoid finite-size effects in the calculations. We show this by 
performing exact diagonalization of different clusters with up to 16 sites. The calculated loss function 
for the extended Hubbard model shows good agreement with experimental spectra. Furthermore, a 
qualitative description of the optical conductivity is obtained with the same Hamiltonian, and the 
same set of model parameters. The comparison with the experiment shows that interladder hopping 
is of minor importance for a realistic description of charge excitations in NaV20s. We find that the 
character of the excitations depends strongly on the direction of momentum transfer. 

PACS numbers: 71.27.+a, 71.45. Gm, 71.10.Fd 



I. INTRODUCTION 

The insulating system a'-NaX^Os belongs to the fas- 
cinating class of highly correlated low-dimensional elec- 
tronic systems. Recently, its physical properties have 
been intensively investigated theoretically as well as ex- 
perimentally. The magnetic susceptibility of NaV205, 
can be well described by a S = 1/2 antiferromagnetic 
Heisenberg chain with exchange interactions oLJ = 440K 
and 560K for temperatures belowtH and aboveB the tran- 
sition temperature Tc ~ 34K. The opening of a spin gap 
at the phase transition is accompanied by unit cell dou- 
bling in the a and b direction, and quadrupling in the c 
direction.13 Based on an early X-ray studyu which pos- 
tulated two different V sites, i.e. magnetic V 4+ chains 
along the b axis separated by non-magnetic V 5+ chains, 
a'-NaV205 was initially identified as an inorganic spin- 
Peierls (SP) materiala similar to CuGeC>3. However, ac- 
cording to a recent crystal structure analysis,!] at room 
temperature all V sites are equivalent with a formal va- 
lence (4-4.5. Therefore, the V sites form a quarter-filled 
laddem which can be mapped onto a S = 1/2 chain,Q 
and the resulting exchange couplings J along the chain 
agree well with the experimental data. In principle, this 
effective chain may show all features of an ordinary one- 
dimensional S — 1/2 Heisenberg chain, including the 
spin-Pcicrls transition. On the other hand, various ex- 
perimental results for NaX^Os indicate a more compli- 
cated transition at T c : The BCS ratio of 2A/kT c for the 
dimerization gap A at zero temperature has the anoma- 
lously high value of 6.44. The. phase transition consists 
of two very close transitions^ and two inequivalent V 
sites were detected by NMRp and attributed to V 4+ and 
V 5+ . This has motivated both theoretical and exper- 
imental wptk,, proposing a pure charge orderiBg|-(,CO) 
instabilityBO or CO coupled to SP distortion.BEJ Re- 



cently, the space group Fmm2 for the low-temperature 
crystal structure was determined from x-ray scattering 
datajifl and three inequivalent V sites in £bp_.V 4+ , V 5+ , 
and V 4 5+ oxidation states were proposed.Ejiij Since this 
space group is incompatible with the NMR measurement 
the low-temperature crystal structure and-the charge or- 
dering pattern are still under discussionOEa 



Charge ordering can 



viewed as a Wigner- 



crystallization on a lattice,E3 which is caused by long 
range Coulomb interaction. Thus different models iiwith. 
nearest-neighbar Coulomb repulsions wejee^udieollJEJ'EJ 
and an in-lindl3 and a zig-zag structuretlJO for NaV205 
were proposed. The calculation of excitation spectra is 
a powerful test for these models, as they can be directly 
compared with experimental data. 

Recently, we studiedE3 the dynamic dielectric response 
of NaV2 05 using an effective quarter- filled ladder model 
(see Sec. [nj). In the present paper we improve on these 
calculations by treating the complete extended Hubbard 
model. Furthermore, we study the role of the interladder 
hopping in more detail. We analyze possible finite-size 
effects in cluster calculations for NaA^Os, and in par- 
ticular, we discuss the optical conductivity. While in 
previous workEj'Ea different models were used for the cal- 
culation of the optical conductivity and the dynamical 
density correlation function, we describe both with the 
same Hamiltonian, and the same set of model parame- 
ters. 

The paper is organized as follows. In Sec. ||, general 
aspects of our approach are presented. The value of the 
interladder hopping is discussed in Sec. fill Finite-size 



effects of cluster calculations are investigated in Sec. IV. 
Section [V| contains the comparison of our results of the 
dynamic dielectric response with experimental spectra of 
NaV205, and in Sec. VI we discuss the optical conduc- 



tivity. Finally, the conclusions are presented in Sec. VII 
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II. GENERAL ASPECTS OF THE 
CALCULATIONS 

We study the dynamic dielectric response and the opti- 
cal conductivity assuming that the electrons in NaV2C>5 
can be described by a quarter-filled extended Hubbard 
model 

H = - J2 *ij (4,<r c j,o + H.c.J +U^2n n n il + 

(id) 

for the 2D system shown in Fig. |l|. (i, j) denotes summa- 
tion over all pairs of nearest neighbors, and spin a =|, j. 



c] are electron creation operators, 



is the occupation-number operator, and [/ denotes the 
Coulomb repulsion between electrons on the same site. 
The hopping parameters ty, and intersite Coulomb in- 
teractions Vij are defined in Fig. |l]. Using second-order 
perturbation theory, one can transform Hamiltonian (|l|) 
into an effective t-J-V model 



n = - ^3 ( e l<A* + h - c - 



(id)," 



+ Yl Jl >i ( Sl ' S J ~ 2 UiU J 1 ) + XI r ' ( 2 ) 



(id) 



(id) 



where c\ a = cJ CT (l — n^-a) are constrained electron 
creation operators, and Si denotes the spin-i operator 
at site i. The exchange interactions between vanadium 
neighbors are parameterized as Jy = 4t?-/J7. 




The loss function in EELS experiments is directly 
proportional to the dynamic density-density correlation 
functionHj By including the long-range Coulomb interac- 
tion in the model within a random-phase approximation 
(RPA) one finds for the loss function 



L(u, q) = Im 



-1 



where 



i f°° 



« q X°(w,q) 



(3) 



(4) 



is the response function at zero temperature for the short- 
range interaction models (|l|) and (Q). depends on the 
energy loss u) and the momentum transfer q. |0) is the 
ground state, p q denotes the Fourier transform of rii, and 
t> q = e 2 A r /(eoe r wq 2 ) is the long-range Coulomb interac- 
tion with unit cell volume v. N is the number of electrons 
per unit cell, and e is the permittivity. The real part e r 
of the dielectric function can be obtained from the exper- 
iment. In the case of NaV2C>5, one finds for momentum 
transfer in a direction e r = 7, and in b direction e r = 5.l3 
The response function Xp( w : Q) an d the optical conduc- 
tivity a a {ui) are connected by the Kubo-Green relation 



a a {u)) = e uj limlm [«qX°(w,q)] 



where at zero temperature a a {ui) is defined by 



(5) 



1 f°° 

fo(w) = —Re / dte lut {Q\j a {t)j a \Q). (6) 
w Jo 

Here j a with a = a,b are the components of the current 
operator parallel to a or b direction. 

Equations (^), for the response function and 
the optical conductivity are valid for zero temperature, 
whereas the experiments have been carried out at finite 
temperatures 




sver, for NaV20s— one finds both 
experimentallyB3'E3'E3 and theoreticallyej that the spec- 
tra depend only weakly on temperature. Therefore, we 
may restrict ourselves to zero temperature. Equations (|^) 
and (^) are evaluated by direct diagonalization using the 
standard Lanczos algorithmEJ which is limited to small 
clusters. For that reason, one can only observe localized 
excitations, and the calculation of the loss function (j^) 
is limited to momentum transfers q > 2it/L, where L is 
the cluster size in q direction. 



III. LARGE OR SMALL INTERLADDER 
HOPPING 



FIG. 1. Schematic structure of the (a,b) planes in NaV^Os 
where the circles denote the vanadium sites. The black sites 
define clusters used in the calculation. 



LDA band structure .-calculations^] and estimations 
based on empirical rulesQ lead to similar values of the 
hopping amplitudes t a , if, perpendicular to and along the 
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ladders. However, there are significantly different esti- 
mations of the interladder hopping t xy . A small[-,value 
of t xy — 0.012eV is found in LDA calculationsfl Ad- 
ditional evidence for a small t xy follows from the weak 
ion along the a-axis as observed by neu- 
On the other hand, a much larger value 



with the above values of and U from Ref. 



The 



magnetic disper. 
tron scattering! 

of t xy — 0.3eV was found from an estimation in Ref. |t]. 



A relatively large interladder hopping of t x 



0.15eV 



was also obtained from a comparison of calculated and 
experimental optical conductivity.L2 ._. 

In Ref. ^9|, a finite temperature Lanczos methodLil has 
been used to calculate the optical conductivity. The 
Lanczos approach is limited to small systems, and a clus- 
ter consisting of two ladders with four rungs [cluster (1) 
in Fig. has been used in Ref. 29 However, when a cal- 
culation is restricted to finite systems, finite-size effects 
may distort the results. Finite-size effects affect delo- 
cahzed excitations, which are important if the hopping 
parameters ty are large, and/or the momentum transfer 
q is small. For that reason, the loss function L(u, q) for 
small q is sensitive to finite-size effects. The calculation 
of the optical conductivity a a (uj) is even more problem- 
atic, as it is obtained from Xp(^,q) in the limit q — > 
[see Eq. (||)]. Since the importance of delocahzed exci- 
tations can be decreased by increasing q, we first study 
finite-size effects of the loss function L(u,q) (Sec. IV). 
Already at finite q we observe large effects even when 
a rather small value of the interladder hopping t xy is 
used. This also means that the optical conductivity for 
the same cluster with a larger value of t xy is even more 
affected by finite-size effects (see Sec. ||). Consequently, 
we can show that the results for the optical conductivity 
from Ref. |2^ are not converged and cannot be used to 
support a large value of t xy . 

In the following, we shall use the hopping parameters 
(t a = 0.38eV, t b = 0.17cV, t xy = 0.012eV) and the on- 
site Hubbard interaction (U = 2.8eV) from Ref. [|. Up 
to now the values of the intersite Coulomb interactions 
V a , Vb and V xy are not known exactly. Therefore we 
choose moderate values for V a and Vb, so that the system 
is close tjO-jLhe-nuantum critical point caused by charge 
ordering.EM3 



IV. FINITE-SIZE EFFECTS 



To investigate finite-size effects, one has to use clus- 
ters of different size. However, the cluster consisting of 
two ladders with four rungs [(1) in Fig. |J is the largest 
to which the standard Lanczos algorithm can be applied 
at present. One way to overcome this problem is to en- 
large the cluster in the direction of the momentum trans- 
fer, and to reduce it in perpendicular direction (so that 
the number of sites is kept constant). Of course, first 
one has to check that one does not distort the results 
by reducing the cluster perpendicular to the momentum 
transfer. To show this we use the effective t-J-V model 



intersite Coulomb interactions are chosen V a — 0.8eV, 
Vb = 0.6eV, and V xy — 0.9eV. First we compare the 
loss function of the original cluster [indicated by (1) in 
Fig. ^, with periodic boundary conditions] with those of 
two smaller clusters (not shown). The first consists of 
two ladders with two rungs for momentum transfer q 
parallel to a direction. The second cluster is an isolated 
ladder with four rungs for q parallel to b direction. Note 
that periodic boundary conditions for these two smaller 
clusters make only sense in direction of the momentum 
transfer. Therefore we choose open boundary conditions 
perpendicular to the momentum transfer. Moreover, for 
the first of the two smaller clusters (momentum trans- 
fer in a direction) one has to use renormalized intersite 
Coulomb interactions V a = V a + Vb and Vb = 2V&. This 
follows from a straightforward analysis of the influence 
of adjacent rungs on the same ladder. 
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FIG. 2. Comparison of the loss function for two smaller 
clusters (full lines, see the text) and the system consisting 
of two ladders with four rungs [dashed lines, cluster (1) in 
Fig. 0| . We use the hopping parameters and the on-site Hub- 
bard interaction of the t-J-V model from Ref. |(| the inter- 
site Coulomb interactions are V a ~ 0.8eV, Vb = 0.6eV, and 
V xy = 0.9eV. The momentum transfer q = 0.3A -1 is parallel 
to a direction (left panel) or b direction (right panel). The 
theoretical line spectra are broadened with Gaussian function 
of width O.leV. 



The results for the different clusters are shown in 
Fig. 2. The loss function for the original cluster (1) of 
Fig. 1 is drawn with dashed lines. The results for the 
smaller clusters with q in a or b direction are shown as 
solid lines in Fig. || (left or right panel, respectively) . As 
can be seen from Fig. [| apart from a small splitting of 
the second peak in the spectrum for q parallel to a di- 
rection there is good agreement between the results for 
the large cluster and both smaller systems. (Note that 
we use a small interladder hopping t xy .) 

Now we can study finite-size effects by increasing the 
two smaller clusters in the direction of momentum trans- 
fer. In this way we obtain clusters (2) and (3) in Fig. [|. 
In the left panel of Fig. || the loss function for cluster 
(2) with four ladders (and for the previous cluster with 
two ladders) are shown for different q values in a direc- 
tion. In the right panel of Fig. || we compare the ladder 
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systems with eight rungs [cluster (3) in Fig. |J and with 
four rungs. One observes large finite-size effects for the 
loss function, especially for small momentum transfer q. 
Note that in particular, the larger clusters (2) and (3) 
from Fig. |l| (solid lines in Fig. |^) lead to additional peaks 
in the spectra. This clearly demonstrates that the origi- 
nal cluster (1) in Fig. [j], which is formed of two ladders 
with four rungs, is not large enough to allow reliable con- 
clusions. 
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FIG. 3. Finite-size effects of the loss function for clusters 
with eight rungs (full lines) and four rungs (dashed lines); 
parameters as for Fig. ^. 
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FIG. 4. Loss function for open boundary conditions (full 
lines) and periodic boundary conditions in q direction (dashed 
lines) ; parameters as for Fig. || 



To test if there are still finite size-effects for clusters (2) 
and (3) , we next compare the spectra for closed and open 
boundary conditions in direction of momentum transfer. 
In the case of open boundary conditions one has to make 
sure that electrons on the edges of the cluster are still em- 
bedded in the local Coulomb potential that results from 
a zig-zag ordered state. For this purpose, sites on the 
edge of cluster (2) in Fig. [j] are assigned an additional 
on-site energy V xy . Sites on the edges of cluster (3) that 
are not occupied in a zig-zag charge ordered state need an 
additional on-site energy Ff,. As can be seen from Fig. |] 
(left panel) there are only small differences between pe- 
riodic (in q direction) and open boundary conditions for 
momentum transfer parallel to a direction, whereas one 
observes systematic discrepancies for small momentum 
transfer parallel to b direction [right panel of Fig. |J. For 
q < 0.2A -1 there are large differences in both direc- 
tions (not shown). We conclude that the convergence of 
the loss function is satisfactory for momentum transfer 
q > 0.2A -1 . In the next section we use clusters (2) and 
(3) of Fig. |l| with periodic boundary conditions parallel 
to q direction for the calculation of the electron-energy 
loss spectrum. 



V. ELECTRON-ENERGY LOSS SPECTRUM 



For the comparison with the experimental loss 
function^ of NaV 2 5 we use the full t-U-V model (@). 
The hopping parameters and on-site Hubbard interac- 
tions are taken from Ref. ^. The values of the inter- 
site Coulomb interactions V a = 0.8eV, V& = 0.6eV and 
V xy = 0.9eV have been adjusted to obtain correct peak 
positions of the loss function. 

In Fig. H the obtained loss function for the full model 
(^) is compared to the experimental spectra taken from 
Ref. |22|. Very good agreement for the loss function with 
momentum transfer parallel to a direction [see panel (a) 
and (b) of Fig. |) is found. In particular, the increasing 
width of the experimentally observed structure at l-2eV 
with increasing q is reproduced. On the other hand, for 
momentum transfer parallel to b direction [panel (d) of 
Fig. H the observed structure at 1.0-1. 7eV is not broad 
enough compared to the experiment [panel (c) of Fig. [| . 
However, the agreement for momentum transfer paral- 
lel to b direction is significantly better fot,thc full t-U-V 
model (0) than for the t-J-V model @E5 Note that for 
momentum transfer in a direction the results for both 
models differ only slightly. This implies that there are 
different basic mechanisms for excitations with momen- 
tum transfer in a or b direction: doubly occupied sites 
are only important for momentum transfer parallel to b 
direction. 

Now let us discuss the character of these excitations 
in more detail. In the ground state all rungs of the lad- 
ders are found to be essentially singly occupied. Shap. 
the system is close to the charge order transitionj23'E2l 
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one can illustrate the nature of the excitations using a 
zig-zag ordered ground state. 




1 2 3 1 2 3 



Energy loss (eV) Energy loss (eV) 

FIG. 5. Comparison of experimental data for NaV^Os 
(left), taken from Ref. ^2], and the calculated loss function 
plotted with an energy resolution of 0.3eV (right). The hop- 
ping parameters and the on-site Hubbard interaction of the 
t-U-V model are used from Ref. [], the intersite Coulomb in- 
teractions are V a = 0.8eV, V& = 0.6eV, and V xy = 0.9eV. 

Disturbance of the charge ordering by electron hopping 
on a rung [process (al) in Fig. ^| is found to be the ba- 
sic mechanism for excitations with momentum transfer 
parallel to a direction. The excitation energy is dom- 
inated by Vj,, and other rungs of the ladders are only 
involved via VJ, and V xy . Note that process (al) can also 
be interpreted as a transition from a bonding to an anti- 
bonding state of a singly occupied rung.Lj In addition to 
this one electron transition, there are also collective tran- 
sitions that involve two [process (al)+(a2) in Fig. || or 
three [process (al)+(a2)+(a3) in Fig. g] adjacent rungs 
on different ladders. The collective nature of these tran- 
sitions results from the interladder Coulomb interaction 
V xy . The role of V xy can best be illustrated in the limit 
of small electron hopping (see Fig. ^). In this case the 
excitation energies are 214, (4V& — V xy ), and (6Vj, — 2V xy ) 



for the processes (al), (al)+(a2), and (al)+(a2)+(a3). 




FIG. 6. Illustration of the basic excitations with mo- 
mentum transfer parallel to a direction (left) and b direction 
(right). In addition to the local transition on one rung [pro- 
cess (al)], for momentum transfer parallel to a direction one 
also observes collective transitions of two or three adjacent 
rungs on different ladders coupled by the interladder Coulomb 
interaction V xy [process (al) + (a2) or (al) + (a2) + (a3)]. The 
creation of one unoccupied and one doubly occupied rung 
[process (bl)] is the basic mechanism for excitations with mo- 
mentum transfer parallel to b direction. Electron movement 
can lead to an increased distance between unoccupied and 
doubly occupied rung [processes (bl) + (b2), (bl) + (b2) + (b3)]. 

The theoretical spectra for momentum transfer paral- 
lel to a direction [see panel (b) of Fig. ||] mainly con- 
sist of three excitations at 1.2eV, 1.6eV, and 1.8eV that 
contribute to the structure at l-2eV energy loss. With 
increasing momentum transfer q spectral weight shifts 
from the excitation at 1.6eV to the one at 1.8eV. The 
three excitations differ in the contribution of processes 
shown in Fig. ^. The excitation at 1.2eV is rather delo- 
calized and consists mainly of multi-electron transitions 
[processes (al)+(a2) or (al)+(a2)+(a3) in Fig. §. The 
other peaks result from more localized processes that are 
dominated by one electron transitions [process (al) in 
Fig. |(|. In addition, the transition (al)+(a2) in Fig. || 
contributes to the excitation with energy loss 1.6eV. 

The excitations that correspond to the structure at 
1.0-1.7eV for momentum transfer parallel to b direction 
sec panel (d) of Fig. § can be interpreted as transitions 
to states with one unoccupied and one doubly occupied 
rung [process (bl) in Fig. 0]. The excitations differ in the 
distance between unoccupied and doubly occupied rung, 
see processes (bl)+(b2) and (bl)+(b2)+(b3). These ex- 
citations are not influenced by the interladder Coulomb 
interaction V xw since the contribution of the V xy part of 
Hamiltonian (|l|) is not changed by electron movement in 
one ladder of the zig-zag ordered system at quarter filling, 
and an unoccupied or a doubly occupied rung can be eas- 
ily moved along the ladder. Since we can only observe 
rather localized excitations by direct diagonalization of 
small clusters there is a lack of spectral weight in our 
theoretical spectra [panel (d) of Fig.J| compared to the 
experimental data [panel (c) of Fig. pf. 
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VI. OPTICAL CONDUCTIVITY 

Before we discuss our results for the optical conductiv- 
ity, we briefly study their finite-size effects. A good agree- 
ment between calculated and experimental optical con- 
ductivity was, obtained from a finite temperature Lanc- 
zos methodrj In Ref. ^ cluster (1) in Fig. |l| with pe- 
riodic boundary conditions was used, and the parame- 
ters of the t-J-V model (||) were chosen as t a — 0.4eV, 
t b = 0.2eV, t xy = 0.15eV, U = 4eV, V a = V b = 0.8eV, 
and V xy — 0.9eV. However, for this system we obtain 
large finite-size effects as can be seen from a comparison 
of the results for periodic [full lines in panel (a) and (b) 
of Fig. Q| and open boundary conditions [dashed lines 
in panel (a) and (b) of Fig. Q . Consequently, the cluster 
formed of two ladders with four rungs is not large enough 
to obtain reliable results for the optical conductivity. To 
overcome this problem we again use clusters (2) and (3) 



are shown. As a result of finite-size effects discussed in 



of Fig. [I] for the calculation of a a and ct&. (In Sec 
we have shown that a reduction of the cluster perpendic- 
ular to the momentum transfer affects the loss function 
only weakly if the interladder hopping t xy is small. This 
statement remains valid for the calculation of the optical 
conductivity if the direction of the momentum transfer 
is replaced by the direction of the electric field.) 




FIG. 7. Optical conductivity of cluster (1) in Fig. |l] [pan- 
els (a) and (b), parameters of the t-J-V model are used from 
Ref. E9j, cluster (2) and (3) [panels (c) and (d), parameters 
of the t-U-V model as for Fig. |H) plotted with an energy res- 
olution of O.leV for periodic boundary conditions (full lines) 
and open boundary conditions (dashed lines). 



Sec. IV we obtain only a qualitative agreement of the cal- 
culated optical conductivity and the experimental data 
of Ref. 32. In analogy to the loss function, there is a lack 
of spectral weight for <ib, and the ratio r = I a /Ib = 3.3 
of total a and b intensities is larger than the experimen- 
tal value of 2.2. Note, however, that the differences of 
the optical conductivity for periodic (full lines) and open 
boundary conditions (dashed lines) using clusters (2) and 
(3) of Fig. |l| are much smaller than the differences for the 
cluster (1). In contrast to Ref. ^9] we have used a small 
interladder hopping t xy to describe the optical conductiv- 
ity. Consequently, we conclude from the comparison with 
the experiment that interladder hopping is of minor im- 
portance for a realistic description of charge excitations 
in NaV20i_ The same was found in a previous theoret- 
ical studyo of the optical conductivity where a t-J-V 
model was investigated which was slightly modified by 



[yl an additional symmetry-breaking on-site energy. 



VII. CONCLUSIONS 

Summing up, we have calculated the EELS spectrum 
and the optical conductivity for the quarter-filled ladder 
compound a'-NaV20s by exact diagonalization of small 
systems using the standard Lanczos algorithm. Our anal- 
ysis shows that finite-size effects for these calculations 
are very important, in particular for the loss function 
with small momentum transfer and for the optical con- 
ductivity. We minimize these finite-size effects by using 
clusters of different shape depending on the direction of 
momentum transfer and the direction of the electric field 
respectively. The results for the loss function are in good 
agreement with experimental spectra. We also obtain a 
qualitative description of the optical conductivity. The 
comparison with the experiment confirms that a large 
value of the interladder hopping is not needed for a re- 
alistic description of charge excitations in NaV205. We 
find that the basic mechanism for excitations strongly 
depends on the direction of momentum transfer. For 
momentum transfer parallel to a direction three differ- 
ent excitations are observed. They differ in their degree 
of derealization. The collective character of the delo- 
calized processes results from the interladder Coulomb 
interaction. The excitations for momentum transfer par- 
allel b direction can be characterized by the creation of 
an unoccupied and a doubly occupied rung. 
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